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SCMMAKV 

I.uin,<ed volume dvnamie equations are dei ived usinn an enerny- state formu- 
lation. I'his technique re<iuires that kinetic and |K)tential eneifiy stale functions 
be written foi- the physical system bein>; investigated. To account for losses in 
the system, a llaylcif;li dissipation function is also formed. Fsinw; these func- 
tions, a LaRraiififan is foi-med ;uid usinfj l.aiir.anf^e's eciuation, the cH|uations of 
motion for the system are tlerived. Tlie results of tlu* application of this techniiiuc 
to a lumped volume are used to derive a model for the free- piston btirlin^; en- 
fiine. rile model was simplified and projjrammcd on .an analog computer. Itc'- 
sults are Riven compariiiR tlie model response witli experimental data 


IN-ltODFCTION 


•\ Rcneral problem presented to many eiiRineers is- Riven a phj'sical sys- 
tem, model it, simulate it, compare results with the physical system, and 
predict performance. The specific nroblcm ol modeliiiR the systmii usually 
leads to eitiu'r a lumped volume or a distributed parameter model. Once the 
particular model type is selected, the expiations of motion that d'scribt' the 
system must be Rcnerated comliK’ve to the model type 

FiiRineerinR systems are Rcnerally (piite lai Rc. They are usually made up 
of more than one type of system. For example, the StirliiiR FiiRine consists 
of thermotK namic. fluid, and mechanical processes. Thus, the cnRincer must 
not only Rcnerate the cijuations of motion for the diffci ent types of systems 
but also the couplinR tHpiations between the systems. There are many metliods 
useil to derive the ix|uations of motion. For a tumped volume model, partial 
dillerential c(|iiaiions describiiiR the flow process can be taken from reference 1, 
and usiiiR approximations, simplified volume dynamic cHiuntions can b«* rcii- 
erated. This was done for a eompri'ssor in reference 2. For nu-ehanu al sys- 



li'ins f rci'-lK)(ly can In* used. In IhcHc cawca all forces arc summed on 

the Iree Ixulics and the resultinR filiations of motion can be generated In these 
cases however, any combination ol the different ty|H‘s of systems ri>qutres a 
separate generation of the coupling f|uations, A more unified approach is to use 
:ui energy- state lormulation of the system. This formulation (ref 3) requires 
kinetic energy, (Kitential energ> , and loss functions to be written for the whole 
system. These functions are comprised of terms fiom all the various types of 
processes contained in the physical system. Once these functions are formeil, 
a l.agrangian and a Uayleigh dissipation function are formed The I'quations of 
motion lor the system are then derived from these functions by Lagrange's 
fiuation. From this method, the coupling between the various systems falls out 
directly as op|K'a.'d to deriving it separately. 

This papiM- applies (he energj- state formulation to the lumped volume twhnique. 
I'sing the results, the filiations of motion for the free- piston Stirling engine are 
derived. The resultiuit model was then simplified to match an c.xperimental test 
setup run at Lewis Itesearch Center. The model was piogrammed on an analog 
computer. A comparison between the model am', the ex|HM*imcntal data is given. 


Js’IIKLlNC CYCLK 

before describing the energy- state approach to the free- piston volume dynamics, 
a brief discussion of the Stirling cycle will be given. A complete treatment can 
be found in reference 1. Idealr/.ed l’-\’ and T-S diagrams are given in figure 1. 

'fhe cycle consists of two isothermal and two constant volume processes Practi- 
cal implementation ot the cycle uses two pistons separated by a regenerator 
dig 2). A regenerator is a |K>rous material which alternately stores or rejects 
heat to a working tliiid The piston phasing is set by the cycle dynamics :ind the 
relative masses of the pistons. The volume between piston f :uid the regenerator 
IS die expansion volume. This volume is kept hot tlirough heat addition 
on fig 1). The other space is called the compression space and it is kept cold 
f>> heat extraction. The combination of heat addition and extraction maintain a 
linnperature gradient acioss the machine. 

Prom tigure 1, the cycle is generated by 

(1) With piston 2 at lull stn<ke and piston I at minimum stroke, all the working 
lluiti is in tile compression space. The volume is at its maximum value while 
luessiire and temperature ate at minimum values. 

(2> Piston 2 IS then movetl in while piston 1 remains fixed Volume df-reases 
and pressure rises since the working fluid is compressed. TemiH-rature remains 

‘ min ss 1 - 2). 



(H) Tor llu* proot'ss (2 - M) Ui(h pit>ton& arc imnrtl lo thr Iflt in tiurh ;i man 
lu'i- ihal llu' volunu* rrmaint> c’ont»lant. Tlu* ivuriu'i ator nlv*‘s olf lira! lo llio 
oohl working Ihiul, lhu« the UMiipfratim* and prcbbiirf riho. 

(II Wlu-n pislon 2 roaclu't* minimum Hlrokr, it tilopb hut piston 1 i-ontiniU‘s 
to movi'. biiuT vtilumr iiu ivasos, prossiirt' dot-ivasos Invaust* ot ('\panston. 

Hu- ti'iiifuMaUiro ivmains t-onstani sini’o IumI is adiiod (pro<’«'ss .'t - I). 

(5( riiiallv. IkiiIi |iistons movo lourlhor (constant vohimet back to llic original 
|N)sitions. Tilt' hot working lluid ^ivcs ol! heat to the rcucncralor, and prt'ssiirc 
ami Icmpcralui i' lii'crcasc. 

The cycle shown in lin^irc I is ideal. t’icncrall\ the I’veb's an- more I'ounded. 
1 he eyelt' thus accepts lieal, converts some ol the heal lo work, then r»'jeels tlu* 
residual he.-it. 


I’KKK-PIM'ON MllU.lNd KNCINf 

I ii'e piston eniiiiie is one in which the nas ilyiiamies r.itiu’r than mechanical 
linkage determines the relative phase belWiU’n the two pistons - one of whudi is 
called the displacer pislt n anil the other, the jH»wer piston. Much ol the free 
piston Mirlinn enume development is atlribuled lo Ifeale (ref. hi. I'lnure :t shows 
a sehemaiie ol a dual Iree piston engine. The dual as|K'cl ol the engine is to 
prov iih* loree b.danein}' lor llu* lar^e pi>wer pistons. .A nood e\planaiu>n ol llu* 
working ol I iv» piston Mirlinn engines is un eii in releri nee I Itriellv, the 
iMuinei* spaces tuiietion as (ias stiriii^s for the displacer and jiower iwsions. .As 
llu* i*\pansioM space ^:is is healed, pressuri* increases and llu* pistons nunc out 
until their respei’tive iKuinee spai*e pressure lorees (^el laruer th.m llu* eombini d 
driving (iressures and (he piston inerlia lorees The method ol ha\inj; one piston 
move wliili* the other remains slationarv is accomplished b\ hav in>i llu* power 
piston mon ....issivi* ilian llu* dis|ilaei*r piston, riuis, the dis(tlaeer moves much 
laslt*r and movi*s llu* viorkinj; lluid bai'k ,ind lorih throim;h tlu* re|ii*neralor Irom 
tlu* expansion space lo llu* compression spaei* I'lu* res'..lt.inl pressure changes 
in the working space caused by cyclic heal addition ruul rejection together with 
the tuned res(HiMse ol the pressures in llu* bouiue s[viees keep the system running. 

MO Ml ' I 

I'lu eiu'rgv- state lornuilation is apidu d to a luni(H*d \ oluine model ol the 
Iree (iislon .Milling engine. Tlu* formulation reiiuires *hal kim tie, potential, and 
loss i*nergies be derived lor tlu* i nitre system, I'he resullani eiu*rg\ st.ile It'rins 
are used lo lorui a l.agrangian and .i Uavleigh dissipation lunelion. Once ihesi* 
tiuielions are lormed, Lagrange's Ki|uation is used to der ve the equalions ol 
motion lor the system Thesi* energv lerms must be i i *,.*m . ali/td eoorihnaies 
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and Uu> forces aclin({ on the system must be generalized. The generalized coor- 
dinates are those which are mutually inde|H*ndent. The number of generalized 
coordinates is injual to the number of degrees of freedom for the system. Typical 
systems and their respective generalized coordinates are given in reference G. 

For a free- piston Stirling engine, there are thermodynamic, fluid, and mechanical 
pmeesses involved. Their respective coordinates are entropy, volume, and dis- 
placement. 

Thermodynamic processes require the use ol very complex energy- .state 
functions, ileference G gives a presentation of such functions and shows that the 
is'rfect gas law is an incomplete description of the state ol a gas. Rather than 
use the complex energy-state functions some simplifying assumptions are made: 

ll) Since lor this afiplication, the fluid velocities and hence the Mach numbers 
an* low, the compliance of the fluid will be considered as a simple spring, that 
IS, pressure is pro|K)itional to volume. This assumption is justified in reference 7 
for small pressure differentials. 

(21 OiffercntJal tenifierature changes are calculated from a separately im- 
jKised energy e<juatlon. 

VVliile these assumptions seem constraining, it will be shown that the technique 
leads to ;ui equation set etjuivalent to the equations commonly used to model volume 
dynamics. 


btatc Functions 

In order to simulate the fluid dynamics, the flaiil system is considered as a 
collection ol discrete lumps (luminng). t)nce the lumping is done, energy, co- 
energy slate functions and a loss function can be derived in order to use the energy- 
state lormulation. i’hese functions result in inductive, capacitive, and resistive 
terms in the injuations of motion (ref. G). These terms will be derived for a fluid 
dynamic process in the following sections. 

Figure 4 snows a string ot control volumes representing the discretized fluid 
system. Between the lumps a spring is used to represent a capacitive or energy 
storage term and wall friction is used to represent a loss term. Kach flow 
volume has an associateil mass of fluid and a displacement volume How rate 
represented by V and Q resfiectively. 


Kinetic Fnergy 

The kinetic energy for the system results from the fluid inductance. From 
relerence G, the pressure momentum is defined as: 



3 


I’j, “IW (U 

All »ymlx)lt> an* Rivon in Apjx'ntlix A. I ib tlu* fluid iiu'i tin. I'li^un* 5 bhowb Iho 
ktnftic t*iu*ruv and I'oonerpv ficldb. From tlu* fijjurc, tlu* kiiiftic onorny is: 




g'(i” ) d I” 
p p 


when* t'.H* pnnu'b indii-ati* a variabli* of inti’nration. Tlu* kinollc coonerpy lb: 


. 


T,W. = / 
^0 


Pj,(Q') dQ' 


I’binp ('({. (1) and abbuminp a linear inductance: 




Fq. (U pivcb the kinetic coenerpy function needed. Next the inertia mubt 
be licfined. Fipure (> bhowb a plfH.* with flultl llowinp throupli it. I'he fluid is 
driven by the pressure difference. Thus: 

P - p., (5) 

a dl\A,/ ‘ 

The use ol l^/A as the velocity assumes that a well ilcfincd velocitv profile 
is known at each cross section of the pi|u*. .\lso the fluid lumps aic assumed to 
nu)\’c as ripid Ixulies. Tluis. the inertia Is: 


llcncc, tlu* kinetic cnerpy for a typic;il lunq) of lluid volunu* in tlu* coonlinates of 
this studv is: 


^ Ap 2 



G 


Potential Energy 

Tie |)otential energy is not eabily described. From figure 4, If the lumps of 
fluid m:iss move at different rates, (xitenti.al energy can be stored in the control 
volume in the form of an Increase in pressure. For small density changes an 
average density is sufficient, and a linear capacitor can be useti to model the 
compressibility effects. From reference G, for fluid dynamics: 

P = — V (8) 

K 

wlicre K is the (spring) compliance (reciprocal of the spring constant), V the 
volume, .and P the pressure. Since the equations derived in this analysis will 
be used for a gas, i» is of interest to derive the compliance for a gas under gen- 
eral conditions .an', compare the result with cq. (8). 

For a general thermodynamic polytropic process: 

PV" = Constant (9) 


Wfferentlatlng; 


V"d.> + nl’V""^ dV = 0 

dl-.ili’dV 

\’ 


( 10 ) 


This relation is shown in figure '4. Comparing eqs. (8) and (10) gives: 


K 



( 11 ) 


Next, for :in average pressure .and volume, a straight line approximation is 
drawn. Thus: 


K = 


\’ 


n P 


Figure 7 can be inverted and translated as is done in figure 8 to show the fx)tential 
energy and cocnergy fields. For an isotheiTn.al process, n equals 1 .and: 

K = -^ (1»- 

P 


H the process is isentroplc, then n = -y , the ratio of specific heats and: 
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V 1’ 


If displacement coordinates this is equivalent to the reciprocal of the spring 
constant nonnally seen for gas volumes for free piston engines, for example; 

K' - - 

\’ 


This constant is given in reference 4 where K' is the spring constant, hince; 

AQ » K Ap 


and 


A V K A p 


The potential energy using figure 8 is 


-AV 


Vf(AV) = 


AP'(A\”)dAV' 



or 


(18) 


V = J- 

'f ~ 

‘ K 2 


(H) 


Finally, n is taken as 1 when formulating the |X)tential cncrg>- terms for 
the flow piocess. 


Hcsistance 

A Rayleigh dissipation function can be written for Ihe flow losses in the sys- 
tem. In figure 4, the losses are indicated as wall frictions. This behavior can 
be expressed as: 

AP = KQ (!•'>) 

where It is the fluid resistance. The defining fields for the Rayleigh and 
CO- Rayleigh dissipation functions are shown in figure 9. For the Rayleigh func- 
tion: 
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J 0 


rainR eq. (15); 


F 


f 



( 18 ) 


(17) 


ANALYSIS OF STIRLING FRKF PISTON SYSTFM 

The syblem analyzed in this pajK-T conta.'nb fluid, mechanical, and thermal 
elementb. The ,n-evioub bcctionr have determined the varioub energy forms for 
I’rc fluid system. Here these results will be combined with the mechanical 
counterparts and later with the energ>’ equation (thermodynamic elements) to 
yield the complete system eejuations. 


Combined Mechanical and Fluid Systems 


Now, the energy formulation is a|)plied to a combined mechanical and fluid 
dynamic system. Figure 10 shows two pistons exerting forces on a contained fluid, 
for a linear mcchrmical system, the kinetic :uid (xitential energies and the loss func- 
tion are well known and are: 


m 



i=l 




V 


m 



m 



i=l 




(18) 
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Since the displacement coordinate for the fluid system Is volume, here let. 


£ = A X, 
1 Pj I 


( 19 ) 


4 , “ A X, 
i Pj I 


Substituting into ef). (18) yields 


m 


T* 

m 


1=1 


-V • <> 

2 A“ 

J Pi 


m 


■\ *> 


i=l 


2 A“ 

J f^l 


m 


'm = 


i=l 


Pi 
2 A“ 

J *^1 


( 20 ) 


Now for the complete system on figure 10, the total kinetic, ixjtentlal and dis- 
sipative energies for the fluid and piston are given by using eqs. (7), (14), 

(17) and (20) with v^^ = 0 
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Note that in the fxjtential energy function that the spring compliance K has the 
square of the area in it. 


Lagrangian 

The Lagrangian form (ref. .1) for the energy forms used arc given by 
The Lagrangian is 
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LaRranRc's e<|uation ts 


(It 


aoj * * 




specifically 


dc 

1 » 1,2. ,.n +2 

''r v„.4j,42 


( 23 ) 


The system of equations represented bv po o .u 

the »y„om and are the malhcmatlcal deacr’lptlon oMhe n^7"T'' 
havlor. The various partlals In eo ^ ^ ^ mechanical be- 

tc r ms determined by aSVaaj are etermined. The Inductance 


2 
Pi 


^^1 a 2 


a^ _ ^^2^ 2 

al2~"X" 

^=1 

3Q, 


P9 


1^1 


( 24 ) 


if- = 'n«n 


similarly, the capacitive terms, ^ ar 




e given by 
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~ 

34, K 


3V, 


1-1 

^^2 ^ iv -2 


<vi:v /ziLhi 


■R 


1-2 


1-1 




av 


For the resistive term 




K 


i >— 2 


K 


n-l-n 


aF _ ^1^ 


a4 


1 a: 


^F _ ^2^2 

^h' a2 
|)«> 


3F 

aQ , 


= 


— = iry„ 
ay„ " " 


Next, uslnp Lagrange's equation (eq. (23)), for *he 4j coordinate 


ilL , 


Pi 


1 


(25) 


(26) 


( 27 ) 


Finally, substitute in eq. (19) and multiply througli by the piston area 
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MiXj 


(a X, - V.^A 

I Pi 1 V 

I, + -1 — i L L ♦ D 


K 


1-1 


X s 0 

r'l 


since is prcsbure, we can w’lite 
K 




( 28 ) 


E'q. (28) ib the equation of motion for the pibton. Note that the force P.A 

^ Pi 

op|K>bCb the pibton motion and ib the coupling force between the me.'hanlcal and 
fluid bybtenib. Thib coupling comeb naturally from the energy-btnte formulation. 
While thib example ib bomcwhat bimple, it btill illubtratcb the value of the method 
when more than one type of bybtem ib contained in a phybical procebb. 

Finally, if a force ib exerted on left pibton, the equation of motion would be 
the bame except for the generalised force. If a force f ib exerted on the pibton 
in the X direction, then all tlic generalized force componentb are 


qv = 0 




= 0 


qtr = 0 




The generalized forceb for the t>quationb of motion are 

^ d.\, d\- dy 

^ — + q,. ■*■••• — - 

* -'l ao ^2 90 'laOj 'naOj 


(23) 


(30) 


However, bince all generalized coordinatcb are independent, the generalized 
forces are the actual forceb in their coordinate system 


= f/A 


Pi 


^2 = • • • 0 


( 31 ) 


Thus, the equation of motion for the piston wi*h force f in the x direction is 


( 32 ) 



H 


M.x, f I),x, = f - I’, A 
11 11 ^ I’l 


(321 


riu* iH|iiaiion of iiiotlon for tho c<K*iillnato it> 


Mj ^2 * " 2^2 - 


( 32 » 


For Uu‘ flow I'oordinato V 


IfCij . 


IV, - V.,) 14i-Vj) 


S .2 


K 


• llj^i = 0 


1.1 


I tsinn till' avi'i ani’ driibity approximation, imj. (3 1) bi'i’onu'B 

r 


\_ 

AjK 


"l * '*l"l ” ’’l ■ *2 


(31) 


i 35 ) 


I'luib, till' now ixpiationb dorivi'd from tlii' I.a^^rancian formulation an* fully 
I'Muivali’nt to tliobi* wliii'li would bi* dorivi'd from diffi'ri'i.tial nu'tliodb. 


1-llKK I'ldTON Kgi’A I IONti 

A bi’lii'inatii' of tin* fn*i' pibton I'liniiii* ib bliown in fif^uri' li. fiiily half of 
till' I'linitu' wab modi'li'd dui* to rom|Hiiii*nt roubtraintb on tlu* analog i'omputi*r. 
I’liuii thib liguri', tlii'ii* an* bi" on volumob and two pibtoiib. Tlu's tliori* an* 
I'iglil niuatioiib of motion. I'bing tlio ivbultb of iHjb. (32) and (35), tin* iH|uatioiib 
of motion an* 


Afg 




^W2 M1;,W2= 1'„- p,, 
A.,g 1 


f., 

— ^ w.. + U.'.w = I’ - p 
A..g ■’ •’ “i 't-. 


( 


V-", 


(3(i) 


(37) 


(3S) 


( 31 ») 
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^''5*'»k"5-'’l-Wl) 


. "(i *' *Voi.i)“ *V 

*'i;8 


'"‘.l\l ^ ‘\rSl “ ^ »V<''d - ^ch> 

''pS. ^ *’pS.= '’pa/^M - ^,>.-> - - ‘V\, 


(10) 

(41) 

(42) 
(45) 


Tin* I'liorjiy iH|ualion is ik'rivi*il in Ap|H*iulix H. Tlu' nriu*ral form of tho 
iloriviul i>(|uation is 



“in"’ 'm - 


T) 


\> ,(y r . - T) 

out ' out 



WlUMO 


(ID 




upstroam 


for |K>sitivo \v.|^ 
for lU'natiM* \\v^^ 


T 


out 


T 


tiownsi roam 


for |H)sitivo w 


out 


for lU'nali' o w 


out 


TIu' roasoM for tlio switohinn is tiu' Didirootional How oharaotoi istio of tho 
maohino. Thoso soxoii tompor aturo ixiuations oorros|H>nd t«> tlio sovon volumos 
o«>nsi(lorod. Tlu* tomporaturos woro oouploil tt» tl;o How iH|uations tlirounii tho 
itloal Has law. 

Finally. lw>th roHi'iiorator soHinonts woro moiloloil as thoiinal la^s. For 
o\amj)lo lor tho first rogonorator 



. r 

I ^ " . 

^ ' m ‘’in 




( 15 ) 


anil w 


m 


m 


whoro o 


aro do|)oiulont on tin- rogonorator matorial. 
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SlMPI.inm l Uth I’lbTON FNCINK MUDFL 

The Irer piston onjjine shown in fitiuro .'t was run at Fowls llosoarch C’ontor. 

For part of Iho lt*st, the |)ower pistims were Imked in their lorward most positiuris, 
and the displacers were sii’oked up and down by oscillating the pressures In their 
bounce spaces. This was done to study the thermodynamic flow process in the 
machine. 

The simulation was set up to match this ex|K‘riim<nt. This was dom* by using 
only the lem(X‘rature and flow liquations. The displacer was externally driven, thus 
eqs. (42» and (-l.'i) were not programmed. 1 urther, since ihe induciiuice terms 
were very small, they were set to zero in eqs. (;tU) to (Hi. This also had the 
desirable etiect ol reducing th.? number of imalog conqxinentb. 

Figuri* 12 shows a comparison of the model results with the ex|x'rimcntal 
data. Note that ihe comparison was (juite good lor Ihe pressure profile. For the 
temperature profile, however, the amplitude swing ol ihe wave sha^x* is good but 
the transient comparison is not g(X)d. The reason for this is still not completely 
understood Ixit is believed related to gas thei nux-ouples intermittently touching 
the exjKmsion space wall. The exjxTimental stored mass compares well with the 
analytical. The expi'rimental values were calculated by using the pressure and 
temperature profiles shown, calculating tin* volume of the expansion space as a 
function ol the displacer ^x)sltion, and then applying the fx-rfect gas law. 'I , mean 
value of the stored mass is somewhat higher for the simulation, but that is cue to 
the slightly lugher mean pressure lor the simulation ;uul the discrepancy betveen 
the analog and ex|x‘rimental temjxrature profiles. I’hal pressure and mass 
compare wliile temperature do**s not is j^Kissible because ol llu> large UC values 
of tem}H*rature (not shown). 

ItKMAKKS AND CONt'U’SlDNb 

The energy- slate formulation was used to derive the equations of motion for 
a fluid- mechanical- thermal system, Hy appropriate assumptions, the proci'dure 
was shown to lead to commonly used equations for volume dynamics. The (>ro- 
cedure also afforded a unified approach in deriving tlie »*quationb of motion where 
more ihaa one type of process is contained in Ihe physical sysfeni. 

The equations of motion for the free piston blirling engine were derived 
using the energy-stale approach. These equations were simplified and nrogrammcu 
on an analog computer. Uesults from Ihe experimental data ami simulation agree 
reasonably well. 

The structures which were used in this formulation were chosen to give re- 
sults consistent with those which would lx* derived by standard diffi*renlial tt*ch- 
niques. It is to be hoped that with more expt'i ience with the approacli that some 
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ol tlu> ahhiimptioiih m:ulo in thib pn^K'r nui;lit lx* ri'laxod and y<‘t yield tormu which 
can be readily binuilatcd and extend the validity ranjte of the model. 

Finally, the energy equation wax bcparately deriveii for thib paix-r. It wouhi 
i)e worlhwiiile to atti-nipl a complete Lai;ran(;i:ui approach which would handle both 
the fluid mechanical at>(H'clb and tlie thermodynamic abpc‘ctu ab well. 
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II 
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K 

K' 
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I’ 


Al'l’KNOIX A 
SYMBOl.h 

o o 
area, m“; (ft“) 

o ■> 

displacer area, m"; (ft“> 

9 o 

displacer rod area, m“; (ft") 

9 9 

|)o\vcr piston area, m“; (ft“) 

9 9 

lK)wer piston area, m“; (ft“) 

9 9 

|K)\vc.’ piston area, m"; (ft") 

9 9 

|K)\vcr piston rod area, m"; (ft") 

s)K*cific heat at constant pressure, J/tkjj-K); (IUu/(lbm-”H)) 
sjiecific heat at constant volume, J/(kj;-K); (Htu^ (Ibm-'^ID) 
s|K'cific heat of tlie mesh, .J/(kK-K); (IUu/(lbm-^U)) 
piston friction 
cliang:e in 

Ilaylcijjh Dissipation Function 
Cc»- Hayleinh Dissipation Fimction 
(Jene rail zed force 
force, N; (IbO 

fjravity, 10 (kj;- m)/N- sec"), (32.2 (lbm-ft)/(lbf-sec")) 
enthal()y, J^kn; (lUu/lbm) 

o 

rcut'iieralor heal transfer coefficient, J/(sec- m"-K); 
(lUu/(sec-ft"-‘’l{)) 

lluid inertia, (N-sec")/(m^); ((lbf-sec")/(ft^)) 

meclianical equivalent of heal, 1.0 (N-m/J); (778.3 (ft-lbO/Htu) 

5 5 

reciprocal of compliance, m' /N; (ft /IbO 

*3 3 

spring constant, N' m ; (Ibf/ft ) 

Fasranpan 
length, rii; (ft) 
mass, kg; (Ibm) 
polytropic constant 

9 •> 

pressure, N/m"; (Ibl 'lt") 
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|) 

* 1)11 
’’iMlj 

i< 

' i’ll., 

a 

Q 

M 

U 

U' 

b 


r 

r’ 

T 


i 

r 

\ 


\’ 


w 

w 


w 


»\ 

% 

\v 

\v 

X 

(f 


m 


A 

P 


2 ■’ 

iliKplaiMT iKuiiK’i' space pressure. N 'm"; (lbf/ft“) 

<> 2 

pressure momentum, (N-see)/m*’; Ulbl-see)/ft*') 

*> •> 

IKUver piston Ikmuui' space pressure, N m"; (Ibl/ll ) 

/ 2 

l>«»\ver piston Ixiunee sfxiee pressure. N/ m ; (Ibf fl“) 

9 3 

voluim* How rate, m' /see; (It' /sec) 
lieat, J; (lUu) 

liiMt rate, .l/see, (lUu/siv) 

Henerali/.ed force eom|x)nent 

5 5 

resistance. (N-see)/m , ((lbl-see)/ft ) 

o 2 

transformed resistance, (N-see)/(kn-m“); ulbl-see)/(lbm-ft“)) 
ent rop\ . i)^(kn-K); (lltu '(lbm-*’fO) 

Kinetic Knerny I'unetion, N-m, (Ibf-ft) 

Kinetic C\v Kne Function, N-m; (Ibf-ft) 
temperature, K, (‘*11) 
time, see 

internal enerfjy, .I/kg. (lltu/llnn) 

ilisplaeenu'nt volume, m'\ (It'S 

I’otential Knerny Function. N-m: (Ibl-ft) 

work, .1; (Iltu) 

mass, ki;; (Him) 

stored mass, kji, (Him) 

regenerator mass, kg, (Him) 

mass flow rate, kj^see; (Him/see) 

mass How acceleration, kfysee ; (Him/sw ) 

lim*ar displacenn-nt, m; (ft) 

Hem*rali/ed coordinate's 

linear displaeenu'iit coordinate transformation 
ratio of sped lie heats 
clianni' in 

di'iisity. kn/m"'; (Him ft 'S 


variable of int(*^>ralion 
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d/ (U derivalivc with robpoct to time 

i)/0x partial derivative with respect to x 

bubttcripts: 

1-6 vai table desiKnation 

c compression volume 

COLD cooler volume 
d displacer 

e expansion volume 

r Huid 

II heater volume 

1 intermediate volume 

i index 

in into a volume 

m mechanical 

out out of a volume 

p ixjwer piston 

It regenerator 

It^ first regenerator volume 

It,, second rejjcneralor volume 

mesh 


w 
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AI'l’KNIJlX li 
KNE’HCJY KgL'ATION 

The change of bturod energy foi‘ a conlroi volume is cunbidered. A typical 
volume ib bhown on figure 13. Now for low velocity gab flow, kinetic and (Mtential 
energieb can be neglectcni. Then the change in Internal energ.v ib given by: 

d(w L') = dfr- dW + h.„w, - h ,w , 
b ^ in in out out 

Now the following bubbtitutionb are made, for equilibrium proccbb 

I = C..T, h “ c T, dW = - dV 

\ p , 


giving 


(W 


= 


I) HV 


And bince w = w. - w , and v = — E 
b in out ' 

_P_ 

.1 dt 

It is to be noted tliat, bince the Milling engine flow directionb obcillate (througli 
zero) due to the nature of the cycle, bwitcliing logic on tempciaturcb and flows 
lb neccbbary to implement the energy I'tpiation eorrectly in each flow liirection. 


w w. (y T, - T) - w ,tyT . 
b in'*^ in out ' out 


T. 


V 
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